Electric-field control of spin waves at room temperature in multiferroic BiFeO,3 
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To face the challenges lying beyond current CMOS-based technology, new paradigms for informa- 
tion processing are required. Magnonics^ proposes to use spin waves to carry and process informa- 
tion, in analogy with photonics that relies on light waves, with several advantageous features such as 
potential operation in the THz range and excellent coupling to spintronics^. Several magnonic ana- 
log and digital logic devices^ have been proposed, and some demonstrated^. Just as for spintronics, 
a key issue for magnonics is the large power required to control/write information (conventionally 
achieved through magnetic fields applied by strip lines, or by spin transfer from large spin-polarized 
currents). Here we show that in BiFeOs, a room-temperature magnetoelectric material, the spin 
wave frequency (>600 GHz) can be tuned electrically by over 30%, in a non-volatile way and with 
virtually no power dissipation. Theoretical calculations indicate that this effect originates from a 
linear magnetoelectric effect related to spin-orbit coupling induced by the applied electric field. We 
argue that these properties make BiFeOs a promising medium for spin wave generation, conversion 
and control in future magnonics architectures. 



Magnetoelectric multiferroics possess coexisting mag- 
netic and ferroelectric phases, with cross-correlation ef- 
fects between magnetic and electric degrees of freedom 6 . 
As such, they can potentially be used to control spin- 
based properties by electric fields^, with very low asso- 
ciated power dissipation. This feature appears promis- 
ing not only for spintronics, in which information is en- 
coded by the spin-polarization of the electrical current^, 
but also for magnonics that use magnetic excitations 
(spin waves) for information processing 1 . Indeed, just 
as coupling between magnetic and ferroelectric order pa- 
rameters exists in multiferroics, coupled spin and lat- 
tice excitations termed electromagnons have also been 
demonstrated 9 -. Such mixed excitations exist at low tem- 
perature in multiferroic manganites^ and are suspected 
at room temperature in BiFeO^. Although in the for- 
mer systems, electromagnons have been shown to depend 
on magnetic fields^, their predicted sensitivity to electric 
fields, at the heart of their potential applicative interest, 
remains to be demonstrated. 

In this paper, we focus on BiFe03(BFO), a special case 
among multiferroic materials. Indeed, BFO has a very 
high ferroelectric polarization reaching 100 ^iC/cm 2 ^ 
below Tc = 1143 K and becomes an antifcrromagnet 
below T_/y = 643 K. The spins form a cycloid structure 
with a wavelength of Aq = 62 nm and an associated cy- 



cloid wavevector equal to Q = 27r/Acr2r— . The magnetic 
cycloid lies in the (-12-1) plane, formed by the ferro- 
electric polarization || [111] and the cycloid wavevector 
Q || [10 — 1] (Fig. QJl). An electric field can simultane- 
ously induce the flop of the polarization and the cycloidal 
plane^. The electrical switch of antifcrromagnetic- 
ferroelectric domains in BFO has been observe d ' and 
controlled in a ferromagnet-multiferroic heterostructurei 
Here we report direct electric-field control of spin wave 
states in BFO at room temperature. 

In BFO, inelastic light scattering reveals two species 
of spin wave excitation modes. These are observed in 
Raman spectra as two series of sharp peaks labelled cy- 
clon ((/)„) and extra cyclon modes (t/j n ) below 70 cm , 
the frequency of the lowest phonon mod o 10 ' 19 ' 20 . The 
cyclon cj) n and extra cyclon ip n modes correspond respec- 
tively to the spin oscillations in and out of the cycloidal 
plane (Fig. []Ji) and their Raman scattering signature are 
reported in Fig. [TJ;. These two species of magnons orig- 
inate from the translational symmetry breaking 2 ^ of the 
cycloidal ground state. The spin wave momentum is no 
more conserved, and can be increased or decreased by 
a multiple of the cycloid wavevector Q. This leads to 
magnon zone folding at the Brillouin zone center, where 
<p n and ip n correspond to magnon wavevector equal to 
k = 7iQ. 
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FIG. 1: Electrical control of spin waves, a, Depicts the 
spin cycloid ground state along with its low energy excita- 
tions when the ferroelectric polarization Pi points along the 
[111] direction. The spin excitations correspond to in plane 
modes <j> with an ellipsoidal shape elongated along y (cyclons) 
and out-of-plane modes ip (extra-cyclons) as ellipse elongated 
along the tangent vector which belongs to the xz plane, b, 
Schematic diagram of the experimental set up used to apply 
an external electric field E to BiFeOs single crystals and to 
probe its spin excitations using Raman spectroscopy. The 
electrodes are deposited on the (010) plane. A positive (neg- 
ative) electric field parallel to the [010] ([0-10]) direction al- 
lows the orientation of the ferroelectric moment Pi (P2) along 
the [111] ([1-11]) direction. Raman scattering is performed 
through the indium tin oxide (ITO) top electrode with light 
polarizations in the (010) plane, c, Raman spectra showing 
the magnon modes 02, i>2, 03, tps and 04 when the applied 
field ranges from to 75 kV/cm. Dashed lines are guides to 
the eye following the shifts in frequency of each mode. 




FIG. 2: Spin wave hysteresis loop. Voltage dependence 
of a, 02 and b, ^2 magnon modes, c, (P-E) hysteresis loop 
at room temperature. The shift of the spin wave frequencies 
follows the polarization loop along 1 (red diamond), 2 (blue 
star), 3 (orange circle) and 4 (green square). 



The experimental set-up used to apply an electric field 
and to detect the optical excitations of spin waves is pre- 
sented in Fig. [Tp. The device consists of bottom and top 
(transparent) electrodes that enable the application of an 
electric field along the [010] and [0 — 10] directions of a 
bulk BFO sample. The effect of the external electric field 
on the spin wave modes is tracked in Fig.[Tp. An external 
electric field applied along the [010] direction induces a 
blue shift of the -0n modes (f/ , 2 and ^3) and a red shift 
of the <f> n modes (02, 03, and 04). The promise of this 
electric-field control depends on our ability to selectively 
tune the spin wave states using the polarization hystere- 
sis cycle and the switching of the polarization vector P. 

Figures [5^,b show the frequency shift of the 02 (cy- 
clon) and tp2 (extra-cyclon) modes as a function of the 
polarization cycle. This figure demonstrates that the spin 
wave modes are directly connected to the polarization 
hysteresis loop shown in Fig.J^F- With an initial polariza- 
tion state Pi || [111], decreasing the applied electric field 
(E || [010]) induces an increase of the 02 mode frequency 
(red diamond, Fig. [5^) following path 1 in the polariza- 
tion cycle. Traveling along path 2 with a negative electric 



field (E || [0-10]), the frequency of the 02 mode (blue star) 
increases up to the flip of P along the [1-11] direction 
(P2) at critical electric field E c ps —35 kV/cmi£. Beyond 
E c , 02 's frequency decreases. The hysteresis cycle can be 
completed along paths 3 and 4 leading to frequency shifts 
that are symmetric to the ones along 1 and 2. The 1^2 
mode is also directly connected to the polarization loop, 
but exhibits an opposite behaviour with respect to the 
02 mode (Fig. [3p). The extra-cyclon mode tp2 presents 
a much stronger frequency shift under the effect of an 
electric field of 50 kV/cm, Au;^ = (+2.1 ± 0.1) cm" 1 . 
Remarkably, the ip2 frequency is seen to increase by over 
5 cm -1 at maximum E field, corresponding to a 30% shift 
of its natural frequency. This experiment clearly demon- 
strates the coupling between the spin waves and the fer- 
roelectric character of the material and that a spin wave 
frequency hysteresis loop can be created and controlled at 
room temperature using the ferroelectric hysteresis loop. 

There exists other methods for electric-field control of 
spin waves. Rado et al. detected a shift of the order 
of 3 x 10~ 5 cm -1 (relative shift of 0.01%) in lithium fer- 
rite under the application of similar voltages^. Recently, 
Vlaminck and Bailleul demonstrated a doppler shift of 
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the order of 10 -3 cm -1 (0.3%) using electric currents^ 3 -. 
The present experiment shows that the coupling between 
magnetic and ferroelectric degrees of freedom provides a 
much more efficient method for electric-field control of 
spin waves, with frequency shifts over 5 cm" 1 (30%), 
several orders of magnitude larger than with previous 
methods. 

The strong spin wave frequency shifts that we observe 
may be interpreted as resulting from a combination of 
magnetoelectric interactions. The component of E along 
P (denoted E z ) cannot play any role because of the com- 
petition with the internal electric field parallel to z pro- 
duced by P. Hence we only consider interactions involv- 
ing E±_ — (E x ,E y ). The Dzyaloshinskii-Moriya (DM) in- 
teraction that couples E to M x L is allowed by the R3c 
crystal symmetry of BFO , and also induces a transition 
from the cycloidal to a homogeneous spin state at large 
E±—. The DM interaction produces spin wave shifts that 
are independent of E at low electric field, at odds with 
the data, and that scale as Ej_ at larger field^. Another 
class of magnetoelectric interactions, the so called flexo- 
electric interaction a 25 ' 27 couple E to gradients of L. For 
example, the flexoelectric interaction in a homogeneous 
magnet leads to an instability towards a spiral state and 
to spin wave shifts linear in Ej_ 27 . However, in our case 
of bulk BFO the magnetic order is a cycloid that is itself 
induced by terms that couple P to gradients of L and 
have the same form as the flexoelectric interactio n 25 ! 28 . 
Therefore, there is a competition between these different 
couplings that leads to a renormalization of the cycloid 
wavevector. Since both uj^ and tu^ are directly propor- 
tional to the cycloid wavevector, a large E± will increase 
all spin wave frequencies by the same amount. Thus, the 
DM and the flexoelectric interactions are not able to re- 
produce the experimental data at lower field but might 
play a role at larger field. 

Here, we have found that the spin wave frequency shifts 
can be interpreted using the Landau free energy model 
based on an additional kind of magnetoelectric interac- 
tions induced by the external electric field. Indeed, the 
R3c symmetry of BFO allows the presence of two linear 
magnetoelectric interactions that couple the Neel order 
parameter L directly with E±, 

Here, we have found that the spin wave frequency shifts 
can be interpreted using the Landau free energy model 
based on an additional kind of magnetoelectric interac- 
tions induced by the external electric field. Indeed, the 
R3c symmetry of BFO allows the presence of two linear 
magnetoelectric interactions that couple the Neel order 
parameter L directly with E±, 

Ft = ^ Kl E ± -L ± L z , (I) 

F2 = \^E X ■ [(L 2 y -Ll)x + 2L x L y y] , (2) 

with K\ and k 2 phenomenological coupling constants. 
Physically, Fi and F 2 model additional magnetic 
anisotropy energies induced by E±. Hence their micro- 
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FIG. 3: Electric field dependence of the measured and 
calculated spin wave frequencies. When the applied elec- 
tric field is zero, the energy of the cj> n and ip n modes are 
E ( j >n = nvoQ and E^ n = y/T+ n 2 voQ respectively, with index 
n labelling the modes, no = 1.4 x 10 6 cm/s the fundamental 
cyclon spin wave velocity, and Q the cycloid wavevectopi^. 
The figure shows the spin wave mode energies normalized to 
voQ as a function of the ratio of the applied electric field E di- 
vided by the electric field E^s required to destroy the cycloid. 
When E > Ehs the spin ground state becomes homogeneous 
in space. The closed symbols are the experimental data ob- 
tained following path 1 (insert) in the polarization cycle using 
E || [010] from to 125 kV/cm : 4> 2 (circle), </f> 3 (diamond), cj> 4 
(down triangle), ip% (star), -03 (up triangle). (P is parallel to 
the [111] direction and no polarization flop occurs. Solid lines 
are explicit numerical calculations obtained from a dynamical 
Ginzburg-Landau theory based on the linear magnetoelectric 
interaction Fi. 



scopic origin is related to spin-orbit coupling. The first 
interaction {F\) produces a red shift for both ip n and <j) n 
modes when E± is increased. Thus, F± cannot explain 
the experimental observations at low E± (Fig. [2^,b). On 
the other hand, the interaction F% induces a blue shift for 
all ip n modes, and a red shift for <f> n modes with n > 2 as 
found experimentally, see Fig. [3] This interaction drives 
the system into a homogeneous state when E±_ becomes 
greater than a homogeneous critical field E\- IS oc 1/k 2 - 
The calculated curves show that cyclon and extra cyclon 
mode frequencies converge to two fixed points at the crit- 
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ical field E^ s - A comparison between theory and exper- 
iment (Fig. [J) gives an estimate of E hs w 160 kV/cm 
(above the maximum electric field of 125 kV/cm that 
could be applied to our sample without damage). Cal- 
culations are in good agreement with the experimental 
data in the low field region (see Fig. [3]), but a clear de- 
viation occurs at high fields. Interestingly, cyclon and 
extra cyclon modes seem to merge for E > 0.5Eh s - This 
effect and the difference between theory and experiment 
at high fields might be due to other magnetoelectric inter- 
actions such as the DM and the flexoelectric interactions 
discussed above. 

The demonstration of electrical control of spin wave 
states represents a significant step towards making new 
spin wave-based technologies. The cyclon spin wave dis- 
persion can be written as lo^ = v^k, and the extra- 
cyclon dispersion as aty, = v^^Jk 2 + Q 2 . For E = 0, 
the spin wave velocities are equal, v$ = — vq = 
1.4 x 10 6 cm/siS. However, under an electric field of 
100 kV/cm the <f> and i/j s P m wave velocities are tuned 
by Av,p — —6 x 10 4 cm/s, and Av^ = +4 x 10 5 cm/s 
respectively. While the speed of electron spin precession 
in a magnetic field is generally fixed by intrinsic proper- 
ties of the materials considered, our work shows that in 
a multiferroic the spin wave speed can be continuously 
adjusted and in a different way depending on the spin 
propagating mode by the applied electric fields. 

We have demonstrated that the frequency of spin 
waves can been tuned electrically by over 30% at room 
temperature in multiferroic BiFeOs. This electric-field 
dependence mimics that of the ferroelectric polarization, 
providing a handle for the non-volatile, low-power con- 
trol of spin waves. Thus, aside from its potential for 
spintronics 29 and photonics^, BiFeC>3 emerges as an ex- 
citing platform for testing novel magnonic device con- 
cepts, which further confirms it as a key multifunctional 
material for beyond-CMOS technology. 

Methode 
Experimental method 

BiFeOs single crystals were grown in air using Bi2 03- 
Fe2C>3 flux in an alumina crucible as detailed in re££. 
and present a single ferroelectric and antiferromagnetic 
domain state. Polarized optical microscope images of the 
samples show that the crystal consists of one single fer- 
roelectric domain 1 ^. Neutron measurements on the same 
samples have shown the presence of one antiferromag- 
netic domain 1 ^. A conducting Indium tin oxide (ITO) 
layer was grown by pulsed laser deposition on the top of 
the (010) sample surface to apply a uniform electric field. 
ITO has no Raman signal in the frequency range of in- 
terest. All the measurements were performed in vacuum 
at room temperature on several samples. 

Raman scattering was performed on the (010) sample 
surface in backscattering geometry using the 647.1 nm 



laser line. Raman scattering is collected by a triple spec- 
trometer Jobin Yvon T64000 equipped with a CCD. The 
spot size is about 100 ^m 2 and the penetration depth is 
less than 10~ 5 cm. Reproducible Raman measurements 
have been performed on several points on the sample sur- 
faces. 

The 4> n and ip n modes (n index labels the modes 
from their lowest to highest energy) were selectively ob- 
served using parallel and crossed polarizations in the 
(010) plane, respectively (see Supplementary Informa- 
tion). The polarizations are defined with respect to the 
projection of the cycloid wavevector Q || [10-1] in the 
(010) plane (the direction of the wavevector Q is not di- 
rectly accessible in the (010) plane). In order to show si- 
multaneously several <fi n and ip n modes, Fig. [TJj presents 
measurements without specific polarization. 

Theoretical method 

We consider a total free energy given by 
F = F + Fx + F 2 , with F Q the usual model free 
energy of multiferroic BF O 21 ' 25 , 

F = ^L 2 + \M 2 + jL±+ C - £ (VL,) 2 

z=x,y,z 

-a P P ■ [L(V • L) + L x (V x L)], (3) 

and Fx, F 2 given by Eqs. ([IJ and ([2j respectively [we 
omit the purely ferroelectric contributions to Eq. ([3]) be- 
cause they play no role in the discussion below]. Here 
A, r, G, c, and ap are phenomenological constants. The 
ground state is determined by the condition 4£ = 0; 
when E±_ = the ground state is a harmonic cycloid 
Lq = Lo[sin (Qx)x + cos (Qx)z] with wavevector Q = — 
and amplitude L% = (-A + cQ 2 )/G. The effect of E± is 
to reduce the ground state wavevector Q and to add addi- 
tional anharmonic contributions to L at odd multiples of 

Q. When E± is equal to a critical field Eh s — °f e , 
the cycloid wavevector Q becomes equal to zero, signaling 
a transition to a homogeneous state. 

The spin wave excitation frequencies are determined 
from the Landau-Lifshitz equations of motion. After 
combining the equations of motion for the order param- 
eters L and M we get 



~6F 








IE ~ 



(4) 



where 7 is the gyromagnetic ratio. We find wave-like so- 
lutions by plugging in L = L + (<5L)e Mt in Eq. @, with 
(SL) = (p(x)D + ip(x)y. The unitary vector D is tan- 
gential to Lq, i.e. it is orthogonal to Lq but lies within 
the cycloid plane. Plugging <f>(x) = J2 n 4>n<3 m ® x , and 
VK^) = i>n^ %n< °^ x m to Eq. ^ leads to a matrix equa- 
tion whose eigenvalues give the spin wave excitations at 
ft = 0. 



5 



We want to underline that the free energy model for 
BFO Eq. (J3J) does not contain any magnetic anisotropy. 
Nevertheless, the present experiment shows that an ex- 
ternal E field induces magnetic anisotropy linear in E, 



see Eq. ([2]). In other words, the electric field activates a 
latent anisotropy in the material, leading to a large effect 
in the case of Eq. ([2j . 
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